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ABSTRACT 



A method of estimating the depth of homogeneous mixing 
in the upper layer of a two-layer ocean is discussed. The 
method applies during the period of general heat loss from 
the ocean surface and is based upon convective processes 
occurring within the homogeneous layer. 

The model used is one developed by N. P. Bulgakov; 
density and salinity distribution with depth are plotted 
from the surface through the halocline. The density scale 
is oriented such that the maximum density for any given 
salinity is read directly. One can thus represent the 
maximum density obtainable (for a given salinity at 
freezing temperatures) by the saHnity plot. By also plotting 
the initial density distribution with depth, and the 
average salinity attained in a homogeneous layer mixed to 
the plotted depth, one can determine the maximum depth 
of convective mixing from the intersection of these curves. 

The model was tested using data taken during the 
winter months over a five year period at Ocean Weather 
Station located in the northeastern Pacific Ocean. 

The conclusions reached show that the model has two 
important features; first, it clearly demonstrates the 
important factors causing convective processes to occur j 
and second, it can be useful in separation of effects of 
the convective process from those of other mixing or non- 
mixing processes in the upper ocean layer. 
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1. Background. 

Knowledge of the vertical thermal structure of the 
oceans is a prime determinant of successful anti-submarine 
warfare. Maximum efficiency in the use of acoustic 
detection devices in the oceans can be achieved only 
through knowledge of its thermal structure. In particular, 
the ability to predict the depth of homogeneous mixing in 
the upper layer would greatly enhance the effectiveness 
of mobile sonar platforms. Toward this goal, investigation 
needs to be made concerning processes which help determine 
the ocean- thermal- structure distribution. In certain 
geographical areas during particular seasons the process 
of convective mixing is predominant in the upper layers of 
the ocean. 

This paper reports tests of a model for estimating 
limits of homogeneous mixing in the ocean due to convective 
mixing during the cooling season of the year. 
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Introduction. 



The mixed-layer depth is defined as the depth to which 
the upper layer of the ocean has been mixed homogeneously. 
There are several causes for vertical mixing in the oceans. 
Basically they can be sub-divided into two categories 

1) mixing due to convective turbulence, and 

2) mixing due to frictional turbulence caused by wave 
action and currents. 

In addition there are two non-mixing types of fluctuations 
of the mixed-layer depth which are associated with 

1) internal waves, and 

2) the vertical convergence -divergence caused by 
changes in the overlying atmospheric pressure 
field, or by convergence-divergence in the surface 
current field. 

The possible influence of internal waves and convergence- 
divergence effects on the mixed-layer depth were not 
evaluated nor accounted for in this study. 

Convective mixing is due to vertical instability. In 
the cooling season there is a net loss of heat from the 
ocean surface caused by upward transfer of sensible and 
latent heat. The subsequent cooling of the water increases 
the density of the surface layer, thereby decreasing the 
vertical stability of the water colximn, or making it 
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unstable. When instability is caused by cooling, convective 
mixing occurs down to a depth where the density is 
equivalent to the increased density of the surface layer. 

This process occurs until there is no longer a net radiative 
loss to the atmosphere or until the freezing temperature 
is reached. 

High rates of sensible heat loss to the atmosphere 
from the ocean surface usually are accompanied by high rates 
of evaporation. Evaporation acts to increase the surface 
density in two ways - by energy loss from the surface through 
the latent energy needed to cause the evaporation; and by 
increasing the surface salinity as the pure water is 
evaporated away, leaving behind in the ocean surface layer 
a greater concentration of dissolved solids. Therefore, 
convective mixing can develop from both sensible heat loss 
and a salinity increase at the surface of the ocean. 

The depth to which convective mixing extends will 
depend to a degree on the initial density gradient present 
in the column of mixed-ocean water. When the initial 
density gradient is small, a given amount of cooling or 
evaporation will set off convection which will reach a 
relatively great depth. When the initial density gradient 
is large, convective mixing will extend downward a 
relatively small amount. For given cooling amounts, a 
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certain critical density gradient can be determined which 
will not permit convective mixing to extend below a certain 
depth. 

The density after mixing of several ocean layers is 
greater than the weighted-average density of the layers 
mixed if a salinity gradient is present, due to the non- 
linear dependence of density on temperature and salinity 
(pressure effects can be ignored when compared to the 
temperature and salinity effects on density in the first 
100m of the ocean below the surface) ClD« This phenomenon 
has been called '’contraction on mixing" by N. N. Zubov [2] . 

N. P. Bulgakov has proposed a model to determine the 
limiting depth of convective mixing for given initial 
conditions this model he has attempted to account 

for this non-linear dependence of density on temperature 
and salinity, especially in brackish water areas in sub- 
polar and polar regions. The model gives a simple numerical 
or graphical determination of the maximum depth of 
convective mixing. According to this model the limiting 
depth of mixing without ice formation occurs when the layer 
has been cooled to its freezing temperature; any further 
increase in convective mixing might be caused by a salinity 
increase due to ice formation at the ocean surface. In the 
greatest part of the open oceans, the upper layer is not 
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cooled to its freezing temperature, and this limiting depth 
of convective mixing usually is not realized. 

Bulgakov’s equation (to be discussed in detail in the 
next section) is the basis for this study xvhich applies the 
convective model to data from a particular ocean region 
during months when convection is anticipated as a dominant 
mixing process. The model gives estimated depths of the 
homogeneous layer which may be verified against measured 
depths. Resolution of consistent differences provides an 
estimate of the validity of the model as well as a suggestion 
of the relative importance of convective mixing in that 
ocean situation. 
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3 . Bulgakov ' s Model . 

In order to understand Bulgakov's approach to 
determination of the limiting depth of convective mixing, 
several definitions are needed. 

1) Density values are in terms of sigma values where 

sigma { O* ) = (density - 1) x 10^. 

2) Reduced density {C^ ) is the maximum density for 
any given value of salinity (at the freezing 
temperature if the salinity is greater than 
24 . 695 or at the temperature of maximum density 
if the salinity is less than 24.695 

3) Convective mixing depth ■ or H) which is the 
depth of the homogeneous layer after mixing due 
to convective processes (less than the depth to 
the bottom of the halocline). 

4) A one-layer ocean is defined as one where the upper 
boiindary of the halocline ( thermocline) reaches 
the surface, with consequent increasing salinity 
(decreasing temperature) commencing at the surface 
and continuing to a finite depth. 

A two-layer ocean is defined as one where the upper 
boundary of the halocline (thermocline) lies below 
a homogeneous mixed upper layer. 
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Zero level (subscript 0) which is the top of the 
halocline in a two-layer ocean; the ocean’s 
surface in a one-layer ocean, 

6) Actual density (cr^) which is the density 

corresponding to the temperature and salinity at 
any given level. 

Bulgakov’s approach relates, for a vertical column of 
ocean extending through the halocline, the profiles of 
initial density distribution and average salinity after 
mixing (plotted against depth) in a basic equation pj . 

Figures 1 and 2 are diagrams of the geometric representations 
for a one-layer and a two-layer ocean respectively. The 
initial salinity distribution is represented by the line 
labeled S; due to the corresponding scale, this line also 
represents the reduced density distribution, C-c , prior to 
mixing. Assuming that salinity varies uniformly with depth, 
the line labeled S represents the mixed-layer salinity 
with changing depth of the convective layer as well as the 
reduced density of the layer after mixing; these values are 
plotted against the existing depth of convective mixing. 

The initial density distribution is represented by the line 
labeled cJt. , The important plots on these diagrams are the 
averaged-salinity curves and the initial density 
distribution curves. The limiting depth to which convective 
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'Figo 2 - Determination of Depth of Convective’ Mixing 
for a Two-layeir Oceans ’ 
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processes can occur is determined by the intersection of 
these two curves. It is at this point that the actual 
density (in the initial density distribution) is equivalent 
to the reduced density (or maximum density) of the 
homogeneous layer after convective mixing. 

The geometrical construction of figures 1 and 2 provide 
a basis for the equation 

“ O"o = + Ss . (1) 

It should be noted that even though (1) is derived from the 
geometrical construction of the figures, each of the 
quantities has a physical interpretation. They are: 

1) (Ju — 0~o is the difference prior to mixing 
between the density at the depth of convective 
mixing, CJia , and the density at the zero level; 

2) Ac is the difference between the reduced density 
at the zero level, O’r^ , and the density at the 

zero level, CJo > (surface or top of the thermocline) ; 
this term represents the temperature reduction 
necessary in order to mix the layer down to its 
maximum depth; 

3) is the difference between the reduced density 
at the zero level, , and the density at the 
depth of convective mixing, . 
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Derivation of the depth of convective mixing, for a 

one-layer ocean based upon equation (1) follows. In figure 
1 it is assumed that the initial salinity distribution 
varies uniformly with depth; therefore, the final salinity 
of the mixed-layer is 



^ 

Equation (1) can be transformed to 






Sh^ s< 



( 2 ) 




'O; 




The quantity is reduced to 



( 3 ) 



Ss = 0.38 4^ He > (4> 

d 2 

where is the salinity gradient in the halocline; 0.3^ 
d z- 

is the change in reduced density per corresponding change 

in salinity. The quantity ^ represents the slope of the 

n c 

dashed curve of figure 1, The constant, 0,3^, is obtained 
by 






1 

0.^8 , 



Equation (1) may now be written as 



dg 

d z. 




^0- + 0.38 4^ He, 

d z ^ 



( 5 ) 



^The constant 0.3^ is determined by Bulgakov; this author, 
using Knudsen's Oceanographic Tables, determined the constant 
to be 0.40 for the salinity range of 32,0 to 34.0 
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from which the limiting depth of convective mixing, , can 
be evaluated if the temperature at the surface (to evaluate 
AG"), and the salinity and the density distribution are 
known. 

The derivation of an equation analogous to equation 
(5) for a two-layer ocean resembles that for a one-layer 
ocean, with a few important differences as figure 2 
illustrates. The derivation of a vertically-averaged 
salinity for the mixed-layer after convection is not as 
simple as in the one-layer case. Since Bulgakov does not 
give a detailed derivation for the two-layer ocean case, 
the following reasonable assumptions are made by this author 
in deriving the final equation stated by Bulgakov. The 
assumptions are 

1) the salinity distribution of the upper mixed layer 
is tiniform with depth, 

2) the salinity varies uniformly with depth in the 
halocline, and 

3) the curve of the average salinity when plotted 
against maximum depth of convective mixing is 
essentially parabolic in shape within the halocline. 

Any weighted-average salinity for the convectively-rnixed 
layer will give a much smaller value for the two-layer 
ocean case than for the one-layer case when convection 
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penetrates to a given depth into the halocline. The low- 
salinity water of the upper layer in the two-layer case 
acts to dilute the convectively-raixed layer. 

It can be shown that the two-layer ocean counterpart 
of equation (4) becomes 



^5 



0.19 



, 

d z. Hj ^ 



( 6 ) 



where Hj is the depth to the bottom of the halocline and 

d3 

is the initial salinity distribution with depth in the 
halocline. The equation comparable to (5) then becomes 

^ H = AC7 + 0.13 li 4* • 

Equation (7) may be solved then for H, which represents the 
maximum change in depth of a homogeneous layer due to 
convective processes when all sensible heat above the 
freezing temperature has been removed from this layer. 

Equation (7) is the key equation in this study of convection 
penetration into the halocline; (7) will be referred to 
hereafter as either Bulgakov’s model or Bulgakov’s equation, 

Bulgakov’s equation may be used as a tool for determining 
maximum changes in the depth to the top of the thermocline 
(assuming that the top of the thermocline corresponds to the 
top of the halocline) , From equation (7) it is seen that 
H, the change in depth due to convection, depends on obtaining 
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1) the density and salinity gradients in the halocline, 

2) the actual and the reduced density at the top of 
the halocline, and 

3) the depth to the bottom of the halocline. 

The salinity gradients can be determined either from 
recent oceanographic stations or from climatology of the 
region. Actual density gradients may be approximated by 
using actual temperatures at the top and bottom of the 
halocline from bathythermograph (BT) observations and 
averaged salinities. Reduced density is evaluated initially 
as a function of salinity (which determines the freezing 
temperature from tables; once determined, the reduced 
density scale may be plotted as an alternate scale for the 
salinity scale as in figures 1 and 2; thus the plot of 
salinity versus depth may be read also as the plot of 
reduced density versus depth. 

Thus equation (7) can be used readily to compute the 
maximum H (due to convection processes) for any area of the 
ocean where vertical salinity and temperature structure is 
known. 



iThe temperature of maximum density was the freezing 
temperature due to the high salinity of the test area. Had 
the salinity been less than 24.695 the temperature of 
maximum density would have been used. 
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4. Application. 

The northeastern Pacific Ocean was the area chosen to 
test the usefulness of Bulgakov’s equation for estimating 
the convective change in the depth of the mixed layer as 
described by the temperature structure { thermocline ) in the 
region. The data were taken by the Canadian weather ships 
which occupy the ocean weather station located at 

50 ° north latitude and 145 ° west longitude. 

The area has two distinct features which make it ideal 
for study. First, the station is located near the center 
of the Alaskan current gyral. Thus the advective change of 
vertical temperature structure due to permanent currents 
can be considered negligible. Secondly, the area is 
characterized by a permanent halocline, located at a depth 
of about 100m at ship ''P' . The halocline separates the 
seasonally-affected upper layer from the nearly constant 
features of the lower layer •’4j. The upper mixed layer has 
an annual average salinity of 32 » 7 average 
salinity at the bottom of the halocline is 33.8%o|53- 

The density structure of the area was considered to be 
a function of the temperature and salinity only. The 
adiabatic changes were considered negligible as were the 
effects due to pressure change with depth, in accordance 
with Sverdrup’s analysis of stability within the upper lOOrn 
of the ocean []J . 
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For a more detailed analysis of the physical structure 
of the northeastern Pacific, the reader is referred to 
Dodimead, Favorite and Hirano[5]. 

The data supplied by the Pacific Oceanographic Group 
of Canada included both oceanographic station data and BT 
observations at ocean station "P" [?3 . The data are 
considered to be extremely reliable. Eight oceanographic 
stations are taken (on the average) each month; and at least 
two BT observations are taken daily at 0200z and 1700z hours. 

The data used were from November, 1957 , through March, 
1962, for the months of November through March. Mean 
salinity gradients for each month were determined using the 
oceanographic stations (See values in Table I). The density 
gradients were determined for each observation using actual 
temperatures at the top and the bottom of the thermocline 
from BT observations and salinities of 32.7 and 
respectively. At the top of the thermocline the reduced 
density is 26.245 gm-cm” 3 . All densities are given as 
sigma values, where 

sigma = (density - 1) x 103. 
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TABLE I - Kean Salinity Gradients 



Month 


dS ( %-9 . \ 


November 


.0144 


December 


.0162 


January 


.0134 


February 


.0266 


March 


.0157 



The Control Data Corporation I 604 Computer v/as used 
along with the observed data mentioned above to solve 
equation (7) for the change in depth to the top of the 
thermocline, H, attributed to convective mixing. The 

density gradient, 4^ , reduced density, cr ^ , and actual 

d z 

density, , were evaluated from the observed data while 

the salinity gradient, 4^ , was taken from the averaged 

d z 

data. After obtaining the quantity, H, this value was then 
added to the depth of the top of the thermocline to obtain 
the predicted depth of the mixed layer. 

As a second application of equation (7), the computer 
was used to solve for the density gradient using arbitrary 
values of H (Im increments) , averaged values for ACT and 

, and arbitrary values for H j . An analysis of various 
terms of equation (7) rewritten as 
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shows that term (1) predominates in determining ^2” , since 

d z 

the vertical salinity gradients are small in the open 

hi 

oceans as is the ratio in term (2) . At ship ’T'" term 

'■’3 

(2) is negligible in comparison with term (1) for the 

computation of ^ , and equation (7a) may be written as 

d z 

^ (7b) 

d-z. H 

Table II is based upon equation (7b) using Im increments of 

H and averaged values of Ad (values are given in the table) . 

This table, therefore, may be used to determine H by entering 

with a computed value of The value of may be 

computed from a current BT observation as outlined previously 

for application with equation (7). The interpretation which 

must be placed on the quantity is that it represents 

dz. 

the critical value of the density gradient; at density 
gradient values greater than , the convective processes 

will terminate at a depth equal to or less than the 
determined value of H. 
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TABLE n- DENSITY GRADIENT VS CONVECTIVE DEPTH CHANGE 



NOVEMBER 


DECEMBER 


JANUARY 


FEBRUARY 


MARCH 




H 


If . 


H 


dz 


H 


£LST 

dz 


H 


ff 


H • 


A O* 


- .961 


ACT* .545 


ACr.= -.509 


ACT. .487 


AO*. 


.489 


. 96 i 


1.0 


.545 


1 .0 


.509 


1.0 


.487 


1.0 


.489 


1.0 


. U 8 0 


2.0 


.272 


2 .C 


.254 .. 


2.0 


.243 


2.0 


.244 


2.0 


.320 


3.0 


. 102 - 


3 .C 


.170 


3.0 


.162 


3.0 


.163 


• 3.0 


.24 0 


4.0 


.136 


4 . C 


.127 


. 4.0 


.122 


• 4.0 


.122 


4.0 


.192 


5.0 


. 109 


5 . C 


.102 


5.0 


.097 


5.0 


,098 


• , 5.0 


.160 


6.0 


.09 1 


6.0 


.085 


6.0 


.081 


6.0 


.081 


> .' 6.0 


.137 


7.0 


.078 


7 .C 


.073 


. 7.0 


.070 


7.0 


.070 


7.0 


.120 


8 . 0 ' 


.068 


8 .C 


.064 


8.0 


.061 


8.0 


.061 


8.0 


.107 


9.0 


.061 


9.0 


.057 


9.0 


. 054 . 


9.0 


.054 


9.0 


.096 


. 10.0 


.054 


10 . C 


.051 


10.0 


.049 


10.0 


.049 


10 . 0 


.087 


11.0 


.050 


11.0 


.046 


1 1.0 


.044 


11.0 


.044 


11 . 0 . 


.080 


12.0 


.045 


12.0 


.042 


12.0 


.041 


12.0 


.04 1 


12.0 


.074 . 


13.0 


.042 


1 3 .C 


.039 


13,0 


.037 


13.0 


.038 >' 


13.0 


.069 


14.0 


.039 


14 . C 


.036 


14.0 


.035 


14.0 


.035 


14.0 


.064 


15.0 


.036 


15.0 


.034 


15.0 


..0 32 


15.0 


.033 • 


15.0 


.060 


16.0 


.034 


16.0 


.032 


16.0 


.030 ■ 


16.0 


.031 


16.0 


.057 


17.0 


.032 


17 . C 


.030 


17.0 


.029 


17.0 


.029 


17 . 0 


.053 


18.0 


.030 


18.0 


.023 


18.0 


.027 


18.0 


. 02-7 


18 . O ' 


.05 1 


19.0 


.029 


19.0 


.027 


19.0 


.026 


19.0 


.026 ■ 


19.0 


.048 


20.0 


.027 


20.0 


.025 


20.0 


.024 


20.0 


■.024 


20.0 


.046 


21.0 


.026 


21.0 


.024 


21.0 


.023 


21.0 


.023 


21.0 


.044 


22.0 


-.025 


22.0 


.023 


22.0 


.022 


22.0 


,022 


22.0 


.04 2 


23.0 


.024 


23 . C 


.022 


23.0 


.021 


23.0 


.021 


23.0 


.040 


24.0 


.025 


24.0 


.02 1 


24.0 


.020 


24.0 


.020 


24,0 


.038 


25.0 


.022 


25 . C 


.020 


25.0 


.019 


25.0 


.02 0 


25.0 


.037 


26.0 


.02 1 


26.0 


.020 


26.0 


.019 


26.0 


.019 


26.0 


.036 


27.0 


.020 


27.0 


.019 


27.0 


.0.18 


27.0 


.018 


27.0 


.034 


28.0 


.019 


28.0 


.018 


28.0 


.017 


28.0 


.017 


28,0 


.033 


29.0 


.019 


29 . C 


.018 


29,0 


.017 


29.0 


.017 


29.0 


.032 


30.0 


.018 


30.0 


.017 


30.0 


.016 


30.0 


.016 


30.0 


.031 


31.0 


.018 


31.0 


.016 


31.0 


.016 


31.0 


.016 


• 31 . 0 , 


.030 


32.0 


.017 


32.0 


.01 6 


32.0 


.015 


32.0 


.015 


32.0 


.029 


33.0 


.017 


33 . C 


.015 


33.0 


.015 


33.0 


.015 


35.0 


.023 


34.0 


.016 


34.0 


. 015 '' 


34.0 


. 014 . 


34.0 


.014 


34.0 


.027 


35.0 


.016 


35.0 


.015 


35.0 


.014 


35.0 


.014 


35.0 


.027 


36.0 


.015 


36.0 


.014 


36.0 


.014 


36.0 


.014 


36.0 


.026 


37.0 


.015 


37.0 


.014 


37.0 


.013 


37.0 


.013 


37.0 


i .025 


> 30.0 


.014 


38. 0 


.013 


38.0 


.013 


38.0 


. 013 - 


38.0 


.025 


39 . O ' 


. 014 . 


39.0 


.013 


39.0 


.012 


39.0 


.013 


39.0 


.024 


40 . 0 ' 


.014 


40.0 


.013 


40.0 


.012 


40.0 


.012 


40.0 


.025 


41 . 0 ' 


.013 


41.0 


.012 


4 1.0 


.012 


41.0 


.012 


41.0 


.023 


42 . 0 ! 


.013 


42.0 


.012 


42.0 


.012 


42.0 


.012 


42.0 


.02 2 


43 . 0 ' 


.013 


43.0 


.012 


43.0 


.011 


43.0 


.01 1 


43.0 


.022 


44 . 0 ' 


.012 


44.0 


.012 


44.0 


.011 


4^.0 


.011 


• 44,0 


.021 


45 . 0 ' 


.012 


45.0 


.01 1 


45.0 


.01 1 


45.0 


.01 1 


45.0 


.021 


46.0 


.012 


46.0 


.01 1 


46.0 


.011 


46.0 


.01 1 


46.0 


.020 


■ 47 . 0 ' 


.012 


47.0 


.01 1 


47.0 


.010 


U 7.0 


.010 


47.0 


.020 


48.0 


.011 


4 8.0 


.01 1 


48.0 


.010 . 


48.0 


.010 


43.0 


.020 . 




. 011 , 


49-0 


..010 


49.0 


.010 


49.0 


.010 


49.0 
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5. Interpretation of Data. 

Interpretation of data is highly subjective due mainly 
to the defining of a thermocline. This author interpreted 
the top of a thermocline as requiring at least a 0.5C 
decrease of temperature in 10m. The bottom was interpreted 
as the location where the temperature decreased less than 
0 . 5 c in 10m. 

A decision had to be made as to when only convective 
mixing is assumed to have occurred. Since the surface 
salinity changes at this station are very small, it can be 
assumed that evaporation effects on density were negligible. 
Convective mixing is assumed, thus, as a function only of 
heat loss in the surface layers during these cooling season 
months. From BT observations, it is extremely difficult to 
read the sea surface temperature to an accuracy greater 
than O.IC. Therefore, it was assumed that when the sea 
surface temperature decreased by O.IC or more in any period 
of time, only convective mixing took place. 
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6. Determination and Evaluation of Results « 

Table III is a condensation of the statistical analysis 
of predicted mixed-layer depth in relation to observed mixed- 
layer depth. Approximately 1,100 BT observations were taken 
during the period of time covered in the analysis. Only 133 
of these observations were used, in the final analysis, as 
representative of purely convective changes, Bulgakov’s 
equation for determining H in a two-layer ocean was used to 
derive the predicted mixed-layer depth. The predicted value 
was arrived at by adding to the mixed-layer depth at the 
time of prediction the computed depth change due to convection. 

Without exception, the arithmetic mean of the predicted 
mixed-layer depth, MLD^j was greater than the mean observed 
mixed-layer depth, KLD. Equation (7) predicts a maximum 
change associated with the removal from the layer of all 
sensible heat above freezing temperature; since this seldom 
occurs, a correction factor was developed to correct this 
over-prediction of convective mixing influence at the top 
of the thermocline. 

The over-prediction is not linear with temperature 
decrease in the surface layer which prevents the use of a 
constant correction factor. However, since the over-prediction 
is related to the retention of sensible heat, a ratio of 
change in sensible heat to initial total sensible heat was 
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TABLE III - Statistical Analysis of Data Results 





November 


December 


January 


February 


March 


^1 


.323 


.666 


.678 


.871 


.926 


^2 


.872 


. 668 


.698 


.876 


.940 


MId 


68 


100 


108 


121 


124 


MLDj^ 


82 


105 


113 


128 


132 


MLD 2 


68 


101 


106 


120 


122 


SD 


11.8 


12.1 


6.2 


7.4 


11.2 


SD^ 


9.1 


14 . S 


7.8 


7.5 


11.2 


SD 2 


9.8 


11.5 


6.2 


7.5 


10.2 


IW 


9.0 


12.0 


12.0 


14.0 


15.0 



^1 



^2 



MLD 

MLD2 

SD 

SD2 

IW 



the linear correlation of predicted MLD to 
observed MLD 

the linear correlation of the corrected prediction 
of MLD to observed MiLD (meters) 

arithmetic mean of observed MLD (meters) 

arithmetic mean of predicted MLD (meters) 

arithmetic mean of corrected RILD (meters) 

standard deviation from MLD (meters) 

standard deviation from MLD^ (meters) 

standard deviation from MLD 2 (meters) 

amp litude of internal waves at the mean observed 
MLD for waves of a period of 7 hours (meters) 
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used to multiply the initially computed value of H. The 
corrected change in depth due to convective mixing is 

H / I - ^ * w> <x.» ~ ~^o ^ 1—1 I I (g) 

Cox^-ecieJ 1= ' Co^npatec! * ^ 

Tjj^ax maximiim annual temperature of the ocean surface 

and Tq is the temperature of the ocean surface at time of 
prediction. A climatological mean of 13. 4C was used 

in the correction factor [i^ . Thus the correction factor 
is a function of Tq only. A correction factor for each Tq 
has been computed and is displayed in Table IV. 

The arithmetic mean of the corrected mixed-layer depth, 
MLD 2 , was in close agreement with the observed mean mixed- 
layer depth, MLD. 

The linear correlation of predicted to observed mixed- 
layer depth both before and after correction are given in 
Table III. With the exception of November, the correlation 
coefficients increased as the cooling season progressed. 

In November, the great difference in the correlation 
coefficients before and after correction appears to involve 
the high temperatures of the surface waters during the 
early part of the month. An analysis of October was dropped 
because heating season effects were predominantly 
superimposed upon the convective or cooling season effects; 
the mixed-layer depth in October was too shallow to escape 
the frictional mixing of wind-induced turbulence. 
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TABLE IV - Initial Temperature vs. Corrective Factor 



Initial 


Correction 


Initial 


Correction 


Temperature , 
To (C) 


Factor 


Temperature 
To (C) 


Factor 


4.1 


.694 


8.2 


. 3^3 


4.2 


.686 


S .3 


.380 


4.3 


.679 


8.4 


.373 


4 • 4 


.672 


g .5 


.366 


4.5 


. 664 


8 . 6 


.35s 


4 • 6 


.657 


S .7 


.351 


4.7 


.649 


8.8 


.343 


4 .S 


. 642 


8.9 


.336 


4.9 


.634 


9.0 


.328 


5.0 


.627 


9.1 


.321 


5.1 


.619 


9.2 


.313 


5.2 


.612 


9.3 


.306 


5.3 


.604 


9.4 


.298 


5.4 


.597 


9.5 


.291 


5.5 


.590 


9.6 


.284 


5.6 


. 5^2 


9.7 


.276 


5.7 


.575 


9 .S 


.269 


5 .S 


.567 


9.9 


.261 


5.9 


.560 


10.0 


.254 


6.0 


.552 


10.1 


.246 


6.1 


.545 


10.2 


.239 


6.2 


.537 


10.3 


.231 


6.3 


.530 


10.4 


.224 


6.4 


.522 


10.5 


.216 


6.5 


.515 


10.6 


.209 


6.6 


.507 


10.7 


.201 


6.7 


.500 


10.8 


.194 


6 . 6 


.492 


10.9 


.187 


6.9 


. 4?^5 


11.0 


.179 


7.0 


.47s 


11.1 


.172 


7.1 


.470 


11.2 


. 164 


7.2 


.463 


11.3 


.157 


7.3 


.455 


11.4 


.149 


7.4 


. 44 ^ 


11.5 


.142 


7.5 


.440 


11.6 


.137 


7.6 


.433 


11.7 


.127 


7.7 


.425 


11.8 


.119 


7 .S 


. 41 ^^ 


11.9 


.111 


7.9 


.410 


12.0 


.104 


8.0 


.403 


12,1 


.097 


8.1 


.396 


12.2 


.089 
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The standard deviations were derived for the computed 
values of MLDj_ and MLD 2 to determine if they are of the 
same order of magnitude as the observed MLD; in Table III, 
it can be seen that the order of magnitudes are related with 
a maximum difference of 2 . 7 m<. 

The computed mixed-layer depth values, both uncorrected 
(IvILDj) and corrected (MLD 2 ), are shown by months and 
compared to observed values of mixed-layer depth (MLD) in 
figs. 3-7» The wide scatter of points in the figures tends 
to decrease as the cooling season progresses. Internal 
waves, it is believed, are the primary cause for the scatter 
of these data points. The amplitudes of the internal waves 
measured in the northeastern Atlantic by Glinskii and 
Boguslavskii range from 8 m to 30m over a period range of 

7.0 hours to 68.0 hours; the amplitude of waves is a 
maximum at a depth range of 100-200m. As an example , the 
amplitudes at 80m for periods of 24 . 8 , 7 » 0 , I 4 . 8 , 12 . 4 , and 

68.0 hours are approximately 8 . 0 , 90 , 9 . 5 , 13 and 22 m 
respectively (]7j . Assuming that these results can be applied 
to a similar region in the Pacific, the standard deviations 
for MLD, MLD^, and MLD 2 in Table III match closely the 
amplitudes of the 7 -hour and the 12 . 4 -hour internal waves 
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PREDICTED MIXED-LAYER DEPTH (METERS) 







Fig. 3 - Observed vs. predictdd-''iViixdd- layer depth for November 
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PREDICTED MIXED-LAYER DEPTH (METERS) 
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• J 

F.lg. 4 - Observed vs. predicted' mixed -layer depth for December 



27 



PREDICTED-MIXED-LAYER DEPTH (METERS) 
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Fig 5.- Observed vs. pr^ddj(;;ted 'mixed-layer depth for January 

2g 



PREDICTED MIXED-LAYER DEPTH (METERS) 
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Fig 6 -Observed vs. predicted mixed- layer depth for February 
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PREDICTED MIXED-LAYER DEPTH (METERS) 

100 110 120 130 140 150 1G0 170 leO 190 200 





1 j I \ \ I I I I f I 
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Fig. 7 - Observed vs. predicted mixed-layer depth for March 
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7. Conclusions and Acknowledgements. 

It is apparent from the data results that Bulgakov's 
model can be useful in determining the maximum or limiting 
depth of convective mixing. 

The main value of his work is that it gives a suggestion, 
in terms of easily-measured parameters, of the mechanics of 
convective activity in the open ocean. Further study might 
be attempted to separate the effects of convective mixing 
and frictional mixing using Bulgakov's equation as a basic 
tool. Present work in the heating season tends to emphasize 
the wind-induced frictional mixing while ignoring the 
increased density (thus convective mixing) effects of mixing 
two or more layers of the oceans. Contraction on mixing of 
several layers should be accounted for in the heating season 
and Bulgakov's equation might be used since contraction is 
automatically accounted for. 

It can be concluded that Bulgakov's equation is useful 
in evaluating processes affecting the thermal structure of 
the oceans. 

I gratefully acknowledge the encouraging help and the 
critical comments of my faculty advisor. Associate 
Professor G. H. Jung. Ky thanks to Dr. Tully of the Pacific 
Oceanographic Group of Canada for the data made available 
for this study. 
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